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We liave performed a Monte Carlo study of the classical XY-model on two-dimensional Sier- 
pinski gaskets of several cluster sizes. From the dependence of the helicity modulus on the 
cluster size we conclude that there is no phase transition in this system at a finite temperature. 
This is in agreement with previous findings for the harmonic approximation to the XY-model 
on Sierpinski gasket and is analogous to the absence of finite temperature phase transition 
for the Ising model on fractals with a finite order of ramification. 
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1. Introduction 

In a series of papers Gefen, Mandelbrot and Aharony examined the possibility 
of phase transitions for discrete-symmetry spin models (e.g. the Ising model) on 
self-similar fractal structures. By applying renormalization-group techniques they 
found that the critical properties of such systems vary with several topological 
characteristics: the fnoninteger) fractal dimensionality D, the order of ramification 
R and lacunarity [ll]. They established that a lower critical dimensionality is not 
defined and found that the Ising model on systems with given D have transition 
temperature Tc =0 if the minimum order of ramification Rmin is finite, and Tc >0 
if Rmin is infinite. The order of ramification R at a point P of a structure is defined 
as the number of bonds which one must cut in order to isolate an arbitrarily large 
bounded set of points connected to P. For two-dimensional Sierpinski gasket (SG) 
Rmin='i and i?m.ax=4 [H- Subsequently, Monceau and Hsiao Q] did a Monte Carlo 
study of the Ising model on fractals with infinite Rmin which had the same fractal 
dimension 1< D <3 but different structure. They found direct evidence of so-called 
weak universality in such systems, i.e. the critical exponents depend not only on 
D, the symmetry of the order parameter and the range of interactions but also on 
topological features of the fractal, in particular on lacunarity which measures the 
deviation from translational symmetry. 

The question remained as to what happens in the case of n-component spin 
models, with n >2, on fractal structures. In [s', '2] a correspondence between pure 
resistor network connecting sites of a given lattice and the low-temperatureprop- 
erties of the spins with continuous symmetry {n >2) on the same lattice |6|] was 
used to conclude that there is no long-range order at any finite temperature if 
the fractal dimension D <2 even in the case of infinite order of ramification. For 
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Figure 1. The first step in creating Sierpinski gaskets. 



two-dimensional SG D = In3/ln2=1.585. Subsequently, in their study of the effect 
of phase fluctuations of the superconducting order parameter in two-dimensional 
Sierpinski gasket wire networks [t^ Vallat, Korshunov and Beck [8] examined the 
XY-model on such a lattice in the harmonic approximation. They found that the 
energy of a vortex is always finite on the two-dimensional Sierpinski gasket imply- 
ing that there is no Berezinskii-Kosterlitz-Thouless transition 0] associated with 
vortex- antivortex unbinding as free vortices are always present. 

Here we present a Monte Carlo study of the XY-model on two-dimensional Sier- 
pinski gaskets described by the Hamiltonian 

H = - J^cosiOi-ej) , (1) 

where 0< 9i <2tt is the angle variable on site i, {i,j) denotes the nearest neigh- 
bors and J >0 is the coupling constant. In the case of XY-model on periodic and 
quasiperiodic lattices the most important excitations at low temperatures are 
the spin-waves and the harmonic approximation to ([I]) used in ^ is quite good 
for the purpose of studying the low-temperature behavior. We compute the heat 
capacity, the helicity modulus and susceptibility for Sierpinski gaskets of several 
sizes. The size dependence of the helicity modulus and its sensitivity to the bound- 
ary conditions clearly indicate that in the thermodynamic limit there is no phase 
transition at finite temperature. 

The rest of the paper is organized as follows. In Section [2] we describe numerical 
procedure used in calculations and present our numerical results. In Section [3] we 
summarize our conclusions. 



2. The numerical procedure and the results 

The procedure which we used to generate two-dimensional Sierpinski gaskets is 
illustrated in Figure 1 , which shows the transition from the zeroth-order Sierpinski 
gasket (the equilateral triangle of side a) to the first order gasket. Our basic unit 
was the fourth order Sierpinski gasket with 123 sites (a gasket of order m has 
A^=3(3'"-|-l)/2 sites) and the list of the nearest neighbors of each site. The higher 
order Sierpinski gaskets were created as subsequent generations of the fourth order 
gasket and we considered gaskets up to order m=7 with 3282 sites. In going from 
the gasket of order m to the gasket of order m+1 the list of the nearest neighbors 
for sites close to the three corners of each of the three mth order gaskets had to be 
modified. 

To study the statistical mechanics of the Hamiltonian ([T]) on Sierpinski gaskets 



we have used Monte Carlo (MC) simulations based on Metropolis algorithm 11 1 



For a gasket of given order the simulation would start at a low temperature with 
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Figure 2. Calculated heat capacity per site as a function of temperature. 



all phases aligned. The first k MC steps per site (sps) were thrown away, followed 
by seven MC links of k MC sps each. The values of k which we used were 120000, 
360000, and 600000, depending on the system size and statistical errors. At each 



temperature the range over which each angle 6i was allowed to vary was ad- 
justed to ensure an MC acceptance rate of about 50%. The errors were calculated 
by breaking up each link into six blocks of k/6 sps, then calculating the average 
values for each of 42 blocks and finally taking the standard deviation a of these 
42 average values as an estimate of the error. The final configuration of the an- 
gles {^1, . . . , 9n} at a given temperature was used as a starting configuration for 
the next higher temperature. Two types of boundary conditions were considered: 
closed, where the three corners of an mth order gasket were considered to be cou- 
pled to each other (then each site has four nearest neighbors) and open, where the 
three corners are uncoupled to each other. 

In Figure 2 we show the results for the heat capacity per site calculated from the 
fluctuation theorem 

_ 1 {H') - {Hf 

where ks is the Boltzmann constant, T is the absolute temperature and (■ ■ ■ ) 
denotes the MC average. The results for the four Sierpinski gaskets were obtained 
with the closed boundary conditions. In the case of heat capacity the results do not 
depend on boundary conditions as illustrated in Figure 3. For comparison, we have 
included in Figure 2 the results obtained for the Hamiltonian ([T|) on three square 
lattices with the periodic boundary conditions. The sizes of the square lattices were 
chosen so that the number of sites is similar to the number of sites for the three 
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Figure 3. The heat capacity per site as a function of temperature for the sixth order SG calculated with 
closed boundary condition (full triangles) and with the open boundary condition (open triangles). 

smaller Sierpinski gaskets. As there is no long-range order in two dimensions the 
heat capacity per site of the square lattices saturates with increasing system size 
(l3l |. In the case of Sierpinski gaskets there is virtually no size dependence of the 
heat capacity per site. This, however, does not mean that the model ([T]) on two- 
dimensional SG also leads to the Berezinskii-Kosterlitz-Thouless (BKT) transition 
resulting from vortex-antivortex unbinding at transition temperature Tc. The peak 
in the heat capacity of the square lattices is above the BKT transition temperature 
where the heat capacity has an unobservable essential singularity [ij]- The peak is 
caused by unbinding of vortex clusters [l^ with increasing temperature above Tc. 
In the case of Sierpinski gaskets the broad peak in C could result from the average 
energy per site {E) changing monotonically from the values near -2J (each site 
has four nearest neighbors) at low temperatures to near zero at high temperatures 
(disordered paramagnetic phase). 

A much better indicator of BKT transition is temperature dependence of the 
helicity modulus. It measures the stiffness of the angles {9i} with respect to a twist 
at the boundary of the system. When the Hamiltonian ([TJ is used to model a set 
of identical superconducting grains coupled via Josephson tunneling, in which case 
9i is the phase of the superconducting order parameter on grain i, the effect of 
applied magnetic field described by the vector potential A is included via shift in 
the phase difference 6i — 9j by 



where the line integral is taken along the line joining sites i and j and $0 = hc/2e 




(3) 
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is the flux quantum [15|. For a uniform vector potential A one has Aij = 27rA 
(rj — rj)/<I>o. Then, the helicity modulus can be defined [3] as 



(4) 

A=0 



where F is the Helmholtz free energy. Using d'^FjdP? = {d'^H/dA^) — 
{{dH/dAf)/{kBT) + {{dH/dA)f /{ksT) one finds an expression for the helicity 
modulus which is analogous to ([2]). For the Sierpinski gaskets A was taken to be 
parallel to one of the sides of the triangles (Figure 1) and for the square latices A 
was taken to be parallel to one of the sides of the squares. Our MC results for 7 ob- 
tained with closed boundary condition for the Sierpinski gaskets and with periodic 
boundary conditions for the square lattices are shown in Figure 4. One expects a 
finite 7 at low temperature when the angles are in an ordered configuration and 



zero stiffness in high temperature paramagnetic phase. Nelson and Kosterlitz 16l | 
predicted a discontinuous jump in 7 at the BKT transition temperature Tc with a 
universal value ^{Tc)/Tc=2/ti, and the straight line in Figure 4 gives the value of 
the jump at different temperatures. For a finite system the jump in 7 is replaced by 
continuous decrease with increasing T which becomes steeper near the (putative) 
Tc with increasing system size, as indicated by the results for the square lattices 
in Figure 4. The two important aspects of our results are: (1) For all lattices con- 
sidered a rapid downturn in 7 starts near the universal 2/7r-line. (2) While the low 
temperature values of 7 for the square lattices do not depend on the system size, 
for the Sierpinski gaskets they clearly decrease with increasing system size and 
the onset of the downturn in 7, which is in the vicinity of the universal 2/7r-line, 
moves to the lower temperatures. This suggests that in the thermodynamic limit 
(i.e. when the order of the SG m —?■ 00) the angle stiffness 7 would be zero at any 
finite temperature implying no BKT transition at a finite temperature. 

We have examined the effect of boundary conditions (closed vs. open boundary 
conditions) and the results for the sixth order SG are shown in Figure 5. While 
the boundary conditions had no effect on the specific heat (see Figure 3) their 
effect on the helicity modulus is dramatic. The open boundary conditions give 
zero helicity modulus at all temperatures within the error bars, which are much 
larger when the three corners of the sixth order SG are not coupled to each other. 
Such a dramatic change in 7 implies that the closed boundary condition introduces 
significant additional correlations in two-dimensional SG fractal structure. Results 
obtained with the open boundary condition reinforce our conclusion about the 
absence of phase transition at any finite T. 

Finally, for completeness, in Figure 6 we show our results for the linear suscep- 
tibility per site computed from 

{M^)-{Ml 

where M is the magnetization of the system when Hamiltonian ([T]) is used to 
describe magnetic systems. For BKT transition Kosterlitz predicted [17] that above 
Tc the susceptibility diverges as x ~ exp[(2 — r/)6(T/rc — 1)""^], with r/=0.25, b 
and z^=0.5, and is infinite below Tc- For finite systems one gets finite peaks in x 
above Tc which increase in height and move to lower temperatures with increasing 
system size, as illustrated by our results for three square lattices in Figure 6. This 
trend is analogous to what one finds for ([T|) in three dimensions where there is 
long-range order below the transition temperature. Therefore the simulation results 
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Figure 4. The helicity modulus 7 as a function of temperature for Sierpinski gaskets and square lattices 
of different sizes. The straight line gives the size of Kosterlitz-Thouless discontinuous jump in 7 at various 
temperatures. 



for susceptibility are not useful in diagnosing BKT transition unless one has the 
results for very large cluster size and attempts to fit them according to prediction 
by Kosterlitz [lOl ]. We note that the susceptibilities for Sierpinski gaskets obtained 
with closed boundary condition are considerably higher than those for the square 
lattices of comparable sizes. Also, the error bars are considerably larger than what 
was obtained for the square lattices. The effect of boundary conditions on linear 
susceptibility for Sierpinski gaskets is shown in Figure 7. The main difference is 
that open boundary condition leads to much larger error bars at low temperatures 
which is consistent with our earlier observation that the closed boundary condition 
introduces additional correlations in the system. 



3. Summary 

Our Monte Carlo results for the specific heat for the XY-model on two-dimensional 
Sierpinski gaskets show that there is no long range order at any finite temperature 
as the specific heat is virtually independent of the system size. Moreover, the 
cluster-size dependence of the low temperature helicity modulus clearly indicates 
that there is no Berezinski-Kosterlitz-Thouless transition in this planar fractal 
structure in the thermodynamic limit. The low temperature values of the helicity 
modulus 7 decrease with increasing system size while they are size-independent for 
square lattices where the BKT transition does take place. Moreover, in the case of 
Sierpinski gaskets the onset of the downturn in helicity modulus near the universal 
2/7r-line moves to the lower temperatures with increasing size of the clusters while 
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Figure 5. The effect of boundary conditions on the temperature dependence of the the hehcity modulus 
for SG of order m=6. 

for the square lattices only the rate of decline in 7 beyond the 2/7r-line increases 
with the system size approaching a discontinuous jump in the thermodynamic limit. 
Our results are consistent with previous findings of Gefen et al. Q who suggested, 
based on their exact renormalization group treatment of the resistor network on 
Sierpinski gasket, that spin systems with continuous symmetry on this structure do 
not display a finite temperaturephase transition. Also, our findings are consistent 
with the results of Vallat et al. [3] obtained for the quadratic approximation to the 
XY- model on two-dimensional Sierpinski gasket, who found that the vortex energy 
is finite and hence free vortices are present at any temperature. 
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Note added in proof 

After this manuscript was accepted for publication it came to our attention that the 
possibility of phase transitions on lattices of effectively nonintegral dimensionality 
was first examined by Dhar (D. Dhar, Lattices of effectively nonintegral dimension- 
ality, J. Math. Phys. 18 (1977), pp. 577-585). In the case of the classical XY- model 
on the truncated tetrahedron lattice (the effective dimensionality 21og2/log5) no 
phase transition at any finite temperature was obtained. 
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Figure 6. The linear susceptibility for Sierpinski gaskets and square lattices. 
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